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Abstract 

We describe an event-based approach to simulate the propagation of an elec- 
tromagnetic plane wave through dielectric media. The basic building block 
is a deterministic learning machine that is able to simulate a plane interface. 
We show that a network of two of such machines can simulate the propaga- 
tion of light through a plane parallel plate. With properly chosen parameters 
this setup can be used as a beam splitter. The modularity of the simulation 
method is illustrated by constructing a Mach-Zehnder interferometer from 
plane parallel plates, the whole system reproducing the results of wave the- 
ory. A generalization of the event-based model of the plane parallel plate is 
also used to simulate a periodically stratified medium. 

Keywords: computer simulation, event-by-event simulation, interference,, 
optics 



1. Introduction 

Maxwell's theory of electrodynamics forms the basis of the understanding 
of the properties of light [1] . The Maxwell equations describe the evolution 
of electromagnetic fields in space and time [1] . They apply to a wide range of 
different physical situations and play an important role in a large number of 
engineering applications. Maxwell's theory describes physical phenomena in 
terms of waves of electromagnetic radiation, yielding a simple explanation for 
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the observation of interference phenomena. For many apphcations, computer 
simulation methods are required to solve Maxwell's equations, the work horse 
being the finite-difference time-domain (FDTD) method [2]. 

In this paper, we present an alternative to the FDTD method. In contrast 
to a wave-based description, our approach uses particles (photons) that in- 
teract with matter. The simulation proceeds event-by-event, that is particle- 
by-particle. There is no direct communication/interaction between different 
particles: Indirect communication takes place via the interaction with matter, 
an interaction that is modeled by means of a deterministic learning machine 
(DLM) [31I11E]- As we show in the paper, our approach is modular and yields 
the same stationary-state results as those obtained from Maxwell's theory. 

In section |2| we introduce the simulation approach and explain how it 
can be used to describe the reflection properties of a single interface, in- 
cluding interference effects (section |3|. The modularity of our approach is 
illustrated by combining two or more DLMs to describe a homogeneous di- 
electric film (section |4]) and a multilayer (section [s]). Finally, we show that 
the basic building blocks can be re-used without modification to simulate 
more complex optical devices such as a Mach-Zehnder interferometer (MZI) 
(section |6]). 

2. Reflection and refraction at an interface 

2.1. Wave theory 

In classical electrodynamics the laws of reflection and refraction are de- 
rived from Maxwell's theory. In case of a plane wave incident on an interface 
between two homogeneous isotropic media with different optical properties, 
there is in general a transmitted wave and a reflected wave. The angle of 
incidence and the refractive indices of both media determine the direction of 
the transmitted and reflected part. Figure [T] shows a schematic picture in 
the plane of incidence. The relation between the angle of the incident ray 6i 
and that of the reflected ray 6r is determined by the law of reflection, 

er = e,. (1) 



The direction of the transmitted wave is determined by the law of refraction. 




Figure 1: Reflection and refraction of light at a plane interface in the plane of incidence. 
The angle of refraction 9t is determined by the angle of incidence 9i and the refractive 
indices of both media rii and n2 (Eq. [2]). The angle of reflection 9^ is equal to the angle 
of incidence 9i. 
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where Of is the angle of the transmitted ray, ni is the refractive index of the 
first and n2 is that of the second medium (see Fig. [T]). 

For lossless, perfectly transparent media, the wave amplitudes are given 
by the Fresnel formulae [1]: 

^ 2ni cos Oi , . . 

" n2 COS^i + TT-l cos6't "' 

T. = '"T"- , ^x, (4) 

rii cos Ui + n2 cos Ut 

n2 cos 9i - ni cos 9t . 

ri\\ = TT^Ih 

n2 cos Ui + rii cos Ut 

riicosdi -rii cos 9t 

R± = ^— ^^±- (6) 

rii cos Ui + n2 cos Ut 

The components of the electric field vector of the incident electromagnetic 
field parallel and perpendicular to the plane of incidence are denoted by A\\ 
and y4_L, respectively. Ty and Tj_ are the amplitudes of the transmitted wave 
and R\\ and R± are the amplitudes of the reflected wave. 
The reflectivity TZ and the transmissivity T are given by 

^ = RF' 

and 

n2 cos 6*4 |Tp 

Equations ([T]) and ^ are both valid for the parallel as well as for the per- 
pendicular part and both components can be treated separately. Combining 
Eqs. (|3]) to ([s]) and simplifying the expressions gives [1] 

_ tan2(6'j - dt) 
^" " tan2(^, + e,)' 

^ (10) 

sm^(6'j + Qt) 

sin 26*4 sin 2^4 

" " sin2(^, + ^,)cos2(0,-0O' 

_ sin 2Bi sin 26*^ , , 

" sin^(^. + ^.)- ^ ^ 
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2.2. Event-based simulation 

We simulate the behavior of reflection and transmission at an interface by 
using an event-by-event, particle-only approach [3llll|5]. Clearly, in an event- 
based model of refraction and reflection at a dielectric, lossless interface, 
there can be no loss of particles: An incident particle must either bounce 
back from or pass through the interface. If such a model is to reproduce the 
results of Maxwell's theory, the boundary conditions on the wave amplitudes 
in Maxwell's theory must translate into a rule that determines how a particle 
bounces back or crosses the interface. In this section, we specify these rules. 

We call an event the arrival of a single photon at the interface. This 
photon carries a message that can be interpreted as phase or time-of-flight 
information. As events occur one at a time only, there is no communication 
between individual photons, but the exchange and the processing of infor- 
mation takes place within the apparatus that describes the interface. An 
incoming photon will either be reflected or transmitted, depending on the 
state of the processing unit. 

For later use, in addition to the input port and two output ports as de- 
picted in Fig. [T} we add an additional input port that captures light incident 
from the opposite direction in such a way that the direction of a refracted 
outgoing particle coincides with the direction of a reflected particle of the 
opposite input port and vice versa (see Fig. [2]). 

The processing unit in this case is a DLM [31 IH E] that takes messages 
from two input ports and sends out messages on either of two possible output 
ports depending on the internal state. The internal state is updated with 
each message that the DLM receives, i.e. it learns from the events that it 
processes. 

The deterministic learning machine consists of three stages (see Fig. |3|: 
The flrst stage receives an input from the n**^ event, in this case the phase 
information of the photon and the angle of polarization w, for practical rea- 
sons encoded as a four-dimensional vector y„ = {{yn)o,\\, {yn)i,\\,{yn)o,±, (l/n)i,±), 
with {yn)o,\\ = cos(0) cos{zu), {yn)i,\\ = sin(0) cos(w), (?/„)o,± = cos(0) sm{zu), 
and {yn)i,± = sin(0) sm{w). Upon arrival of a photon at one of its two input 
ports, the DLM stores the message in its internal register Yj. with k = or 
A; = 1 if the input was on port or 1, respectively. There is also an internal 
vector X = (xq, xi) with Xi G [0, 1], i = {0, 1} and xq+xi = 1. This vector is 
updated for each event received on port k according to 

{Xn+l)i = a{xn)i + (1 - (13) 
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Figure 2: Reflection and refraction of light at a plane interface in the plane of incidence. 
There are two input ports (0 and 1) and two output ports (0 and 1). The angles are 
chosen such that the direction of a refracted outgoing photon coincides with the direction 
of a reflected photon of the opposite input port and vice versa, allowing for interference 
to occur. 
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Figure 3: Schematic diagram of a DLM that performs an event-based simulation of a 
plane interface. There are two input ports and two output ports. The first stage (DLM) 
updates the internal registers Yq, Yi and x according to the input y. The second stage (T) 
processes the information stored in these registers according to a specific rule (Eq. (14)) 
and the third stage (O) prepares the outgoing message and sends it through one of the 
output ports. 



where < a < 1 is a parameter that determines the speed of learning. Since 
{xn)o + {xn)i = 1 for all n, we can interpret {xn)k as (an estimate of) the 
frequency for the occurrence of an event on port k. 

The second stage processes the information stored in the registers Yq, Yi 
and X according to the rule 



/ iyo)o 

(>o)o,± 
(>^o)i,± 
(>^i)o.± 



xl 
xl 
x^ 
x^ 



/ ni(^o)o,!i 

r±{Yo)o,± 
'^±(^o)o,± 



x^ + mYi)o,\\y/x^ \ 
x^ + ti\{Yi)i^lly/x^ 

x^- r||(ri)o,||v/^ 
x^^- r||(Fi)i,j|^/xI 
xJ^ + t_L(Fi)o,±yxI 
x^ + t±(Yi)i^j_y/xl 

xo - r_L(yi)o,±v^ 
x^-r±(Yi)i^^^ ) 



(14) 
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with 




t|| = 
r_L = 
t± = 

Here we have omitted the event label n. 

The third stage of the DLM prepares the messages 

/ t||(Fo)o,||y^- '^iK^Oo.iiv^ \ 



(15) 
(16) 
(17) 
(18) 



(19) 



and 



xJ^ + t||(Fi)o,||v/^ \ 
+ til (Fi) 1,11^x1 



(20) 



ni(^o)o,i 
ni(^o)i,i 

and generates a uniform random number < r < 1. If ||w^ 
final stage sends z„ = w„/||w„|| through port 0. Otherwise it sends = 
w^/||w^|| through port 1. 

From the above construction of the DLM, it is clear that the connection 
between Maxwell's wave description and the event-based, particle-like sim- 



> r, 



the 



ulation model enters through Eqs. (15)-(18), where the expressions in the 
left-hand-sides of the latter are given by Eqs. (|9|-(12). 

2.3. Simulation results 

As a first validation of the simulation model, we simulate a single interface 
with our event-based method. The initial values of the registers Yq, Yi and 
X are chosen randomly, but properly normalized. For each point in Fig. |4| 
100000 events were simulated by sending messages on port with a randomly 
chosen but fixed phase. The parameters were set to a = 0.99, rii = 1, 
and n2 = 1.52. At the end we count, how many events we have detected 
on output port 0, i.e. the fraction of reflected particles. This normalized 
intensity corresponds to the reflectivity TZ. 
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Figure 4: Reflectivity 7?. of a plane interface between homogeneous dielectric media with 
refractive indices ni = 1 and n2 = 1-52 as a function of the angle of incidence 9i. Triangles: 
Polarization angle = 0; Squares: w = 7r/4; Circles: vj — 7r/2. Each marker represents 
a simulation of 100000 events with a = 0.99. The simulation data are in very good 
agreement with the wave theoretical predictions (solid lines). 
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2.4- Discussion 

We have shown that the reflectivity of a single interface can be simulated 
by our event-based approach. The simulation results are in excellent agree- 
ment with the theoretical predictions of Maxwell's wave theory. Features 
such as the polarization dependence or reflection at the Brewster angle are 
faithfully reproduced. However, this good agreement does not show yet that 
our model correctly simulates interference phenomena. Such a demonstration 
is given in section [3] 



3. Interference effects at an interface 

In section |2| we dealt with the case of input on a single input port only. 
Here we consider the case where particles can arrive, one-by-one, on both 
sides of the interface, that is on both input ports. 

3.1. Wave theory 

According to Maxwell's theory, if particles enter on input port with a 
probability of po, carrying a phase (po and on input port 1 with probability 
I — Po and phase 0i the amplitudes on the output ports are given by 



upo 
i(f>i 



(21) 



with _L} and given by Eqs. (15) to (18). For S-polarization, i.e. the 

_L-component {w = vr/2), the normalized intensity on output port or the 
reflectivity 7l± is given by 

,2 sin^ (ei-Ot) sin 2^^ sin , 

sm{9i - 9t)y/sm29i sm29t cos(0o - 0i) r~Pi \ 

-2 . 2fn Vpo(l-Po)- (22) 

sm (fc'j + 9t) 

For P-polarization {w = 0), the corresponding expression reads 



taii^{9i-9t) sm29ism29t 
tan2(^, + ^,)^° + sin2(^, + ^,) cos2(^, - 9^) 



tan(6'i - 6't)v/Sr25~sm2^cos(0o - 0i) r~r, \ 



+ ^ tM0^ + ^t)sM + 9^cos(9,-9^ VPo(l-Po). (23) 
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3.2. Simulation results 



In Fig. |5]we compare the event-based simulation results to the wave theo- 



retical predictions Eqs. (22) and (23). In these simulations, we send, one- by- 
one, photons with phase 0o on port with probability po and photons with 
phase 01 on port 1 with probability 1 — po- Of course, the results depend on 
the incident angle 6i. 

3.3. Discussion 

Our event-based simulation results are in excellent agreement with the 
wave theoretical description. For a single interface and a single input port 
there are no interference effects and the results are in agreement with wave 
theory, independent of the parameter a. However, in the case of two input 
ports where interference can occur, the value of a is important: The wave 
theoretical predictions can only be reproduced by the event-based simulation 
if a is close to one j3l IH [5] . 



4. Light propagation through a homogeneous dielectric film (plane- 
parallel plate) 

A homogeneous dielectric film between two homogeneous media can be 
regarded as two plane parallel interfaces. Figure |6] shows a schematic picture 
of the system. 

4.1. Wave theory 

According to Maxwell's theory, the reflectivity and transmissivity of the 
film are given by [1] 



^ _ ^^2 + ^la + 2^12^23 COS 2/3 
1 + 

'^12'^23 + 12'"23 COS 2/3 ' 

^3 COS tl^tl^ 

rii cos 6*1 1 + r J2'^23 + '^'^12^2^ cos 2/? ' 



and 



with 



na cos ^3 tj^tl^ . 



ni cos 9i - 722 COS 02 

ri2 = ^ ^, 26) 

ni cos Ui + n2 cos 02 



2ni cos 61 

^12 



ni cos 61 + n2 cos^^2 ' 



(27) 
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Figure 5: Single interface between a ni = 1 and a n2 — 1.52 dielectric with input from 
both sides. Single photons arrive on port with probability po, carrying a phase (po and on 
port 1 with probability l—po and phase 4>i. At any time there is at most one single photon 
in the system. Each marker represents a simulation of 100000 events with a = 0.98. The 
phase difference c/jq ~ (f>i is generated from randomly chosen values of 4>o and 0i (which 
are then fixed for 100000 events). Solid circles: Polarization angle -cu — 0; Open circles: 
m = Tr/2. (a): po = 1/2, 9^ = 0; (b): po = 1/2, 0, = tt/I; (c): po = 1/4, 9, = 0; (c): 
Po — 1/4, 9i — tt/A. The simulation results are in excellent agreement with the theoretical 
expressions derived from wave theory (Eqs. (22 1 and (23)), shown here as solid lines. 
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Figure 6: Homogeneous dielectric film built from two plane parallel interfaces. The film 
thickness is denoted by h. For a plate surrounded by air the refractive indices ni = 77.3 = 1 
and 9i — 9^. 

for a S-polarized wave (ccj = 7r/2) and 

n2 cos 61 - ni cos 6*2 

^2 = ^— ^, 28) 

n2 costal + rii cos6'2 

2rai cos , , 

tu = n\ 29 

n2 costal + Hi COS6'2 

for a P-polarized wave (ro = 0), describing the process at the interface from 
medium 1 to medium 2 and analogous expressions for r23 and ^23- The 
variable /3 is determined by 

2tt 

/3 = —722/1 cos ^2, (30) 

with Aq being the wavelength of the incident wave and h denoting the thick- 
ness of the film. 

4-2. Event-based simulation 

We simulate a homogeneous dielectric film by connecting two DLMs, one 
for each interface and each working as described in section [2j For each 
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interface, we use the appropriate expressions for r and t as given by Eqs. ( 15 ) 



(18) and Eqs. d9])-(12), that is we do not use the expressions from Maxwell's 



theory for the film (Eqs. (24)-(30)). We recover the results of Maxwell's 
theory by the event-based simulation without solving the wave equation for 
the film. This modularity allows us to re-use the DLM model of an interface 
for simulating films, multilayers etc. 

We assume that the path of the particles is as depicted in Fig. [7} The 
incident particle is either reflected or refracted as it hits the first interface. 
In case of reflection it leaves the system on the front side of the film. If it is 
refracted, the particle refracted from the first interface acts as the incident 
particle of the second interface. While traveling from the first to the second 



interface it acquires a phase shift exp(i$), with $ = /3 (Eq. (30)), depending 
on the width of the film. Hitting the second interface there are again two 
options; either the particle is refracted and leaves the system on the back 
side of the film or it is reflected towards the first interface again, but this 
time entering the other input port after acquiring another phase shift on the 
way. Subsequently, a refraction leads to an exit on the front side of the film 
and on reflection the particle is sent back to the first input port of the second 
interface. This continues until eventually, the particle leaves on the front or 
the back side of the film. The arrangement of the two DLMs is depicted 
in Fig. |8} This setup corresponds to the behavior described and shown in 
Fig. [7} 

4-3. Simulation results 

With the setup of Fig. [8| we simulate the behavior of homogeneous di- 
electric films. The first analysis considers a plate of thickness h = n2Ao/4 
(quarter-wave plate) under normal incidence for various values of refractive 
indices. The results for 100000 events and a = 0.99, together with the results 
of Maxwell's theory, can be found in Fig. [9| The values of the reflectivity 



predicted by wave theory (Eq. 24 ) are reproduced with high precision by the 
event-based simulation. Figure [lO] shows the reflectivity under normal inci- 
dence of various plates surrounded by air {rii = 723 = 1) depending on the 
thickness of the plate. Again, the agreement with wave theory is very good. 
The variation with the incident angle for different polarizations is shown in 



Fig. [TT] The wave theoretical predictions are all reproduced by the event- 
based simulation, including features like zero reflectivity under the Brewster 
angle for P-polarized {m = 0) light . 
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Figure 7: A homogeneous dielectric film built from two plane parallel interfaces. Upon 
incidence on either of the two interfaces a photon can either be reflected or refracted. This 

can continue until it leaves on any side of the film. Due to the translational invariancc when 
dealing with plane waves, corresponding translated paths can be regarded as superimposed. 
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Figure 8: Arrangement of two DLMs simulating a homogeneous dielectric film. The output 
on port of the second DLM is fed back into the input port 1 of the first machine. The 
symbol <I> denotes a phase shift that a message acquires when travelling that specific path. 



4-4- Application: beam splitter 

Having shown that the event-based model of the homogeneous dielectric 
film works as expected, we now illustrate the modularity of the simulation 



approach by building a 50/50-beam splitter. The setup is shown in Fig. 12 
The incident angles are fixed to 45°, the width h of the plate is set to h = 
Ao/4n2 and the refractive index n2 of the plate is set to n2 = 1.86, such that 



we get a 50/50-beam splitter for -oj = 7r/2 (see Eq. (24)). 

The incident particles on both input ports hit the beam splitter such 
that the direction of the transmitted particle from one port coincides with 
the direction of the reflected particle from the other input port. Due to the 
translational invariance of plane waves, parallel paths can be overlayed, as 
in the case of the plane parallel plate. 

Next, we consider an experiment where we send, one-by-one, particles 
carrying the phase information (pQ {(pi) to input port (1). The probability 
for a message to arrive on port is po and with probability of 1 —po a message 
arrives on port 1. According to wave theory, the amplitudes bo and bi on the 
output ports and 1 are given by [1] 



bo 

h 



t r 
r t 



(31) 
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Figure 9: Simulation results of the reflectivity TZ oi a quarter-wave plate, built from two 
plane parallel interfaces, for normal incidence and for various choices of refractive indices 
rii, 77,2 and 713. Symbols represent simulation results with a — 0.99 and 100000 events per 
data point. Solid squares: ni — = 1; Solid circles: ni — 1 and 713 — 1.45; Triangles: 
Til = 1 and na ~ 1.8; Open squares: ni = rt3 = 1.45; Open circles: ni = ria = 1.8; The 
solid lines are the wave theoretical predictions (Eq. ([24|)). 
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Figure 10: Simulation results of the reflectivity 7?. of a homogeneous dielectric film for 
normal incidence {9i — 0) in air (ni = 77.3 = 1) as a function of its thickness h for various 
values of 712. Sohd circles: n2 = 3; Sohd squares: n2 = 2; Triangles: 712 = 1.7; Open 
circles: n2 = 1; Diamonds: n2 — 1.5; Open squares: n2 = 1.4; Open triangles: 712 = 1.2. 
The parameter a = 0.99 and 100000 events were processed for each set of parameters. 
The simulation results (markers) are in very good agreement with the wave theoretical 
expressions (sohd hnes, Eq. (24)). 
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Figure 11: Simulation data of the reflectivity 7?. of a quarter wave plate as a function 
of the angle of incidence 9i for different polarizations ro. Triangles: Polarization angle 
tn = 0; Squares: tn = 7r/4; Circles: m ~ Tr/2. The parameter a = 0.99 and 100000 events 
were processed for each data point. The simulation results (markers) are in very good 
agreement with the wave theoretical expressions (solid lines, Eq. (24|) 
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Figure 12: Beam splitter built from two plane parallel interfaces. The incident angle is set 
to 45°, the width to h = Ao/4n2, i.e. to an optical path length of a quarter wave and the 
refractive index n2 is set to n2 = 1.86, so that we get a 50/50-beam splitter (m = 71.3 = 1, 
w — 7r/2). There are two input ports and two output ports. Due to the translational 
invariance of plane waves, parallel paths (indicated by dotted lines) can be matched. This 
is also true for paths emerging from multiple reflections within the plate (not shown here, 
compare Fig. [?]). 
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with r and t given by 

(32) 



ri2 + ^^236^*^ 
r = 

and 

^12^236*^ 



t = — ^. (33) 

1 + ri2r23e2^/^ ^ ^ 



The definitions of ri2 and ti2 are given in Eqs. (26) and (27), and ^23 are 
defined analogously. 

Depending on the phase difference 0o ~ 0i the probability for a particle 
to exit on the output port is given by 

1 



l^ol' = ^ + \/po(l - Po) sin(0o - 0i). (34) 

We have run the single-event simulation with a = 0.98 and Nq + Ni = 
100000 events for each value of the phase difference, where A^^o and A^i denote 
the number of events on the output port and 1, respectively. We determined 
the normalized intensity Nq/ (Nq+Ni) detected on output port 0. The results 



for various values of Pq are shown in Fig. 13 The event-based simulation is 



in very good agreement with the behavior predicted by wave theory. 
4-5. Discussion 

We have shown that we can use an event-based simulation approach to 
describe the behavior of a homogeneous dielectric film. With the proper 
choice of parameters we can use this as a beam splitter. This beam splitter 
can be used as building block for optics experiments like the Mach-Zehnder 
interferometer (see section |6]). 

5. Light propagation through a periodic multilayer 

5.1. Wave theory 

A periodic multilayer consists of a succession of homogeneous layers of 
alternating refractive indices, n2 and n^, and thicknesses, h2 and h^, between 



two homogeneous media with refractive indices ni and n4 (see Fig. 14). For 
quarter wave layers (n2/i2 = n^hs = Ao/4) at normal incidence the refiectivity 
for a total of N interfaces is given by fl] 



n 



ni J y-rii 
ni I \ 713 



A 



if N is odd, (35) 
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Figure 13: Simulation results for a 50/50-beam splitter built from two plane interfaces. 
The parameter a — 0.98 and 100000 events have been processed for each data point. 
The phases 00 and were drawn from uniform random distributions. We measured the 
normalized intensity No/{No + Ni), with A^o and A^i being the number of events on the 
corresponding output port. Depending on the probability po of an event entering on input 
port 0, we get different results which agree very well with the wave theoretical predictions 
(solid lines). Triangles: po = 1; Circles: po = 0.5; Squares: po — 0.25. 
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Figure 14: Scheme of a periodic multilayer. It consists of a stack of layers with alternating 
refractive indices, and ns, and thicknesses, /12 and /13 between two media with refractive 
indices ni and 71,4. Drawn here arc only three periods, but the number of periods can be 
chosen arbitrarily. It is also possible to end with the transition from 712 to 714 skipping the 
last layer with refractive index 723. 
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Figure 15: Network of DLMs for the simulation of a periodic multilayer. Here, only a 
single period {ni — 712 — ^ na — >■ n^) with 3 interfaces is shown, but in principle it can 
be generalized to arbitrarily many periods. The first DLM on the left simulates a plane 
interface between media with refractive indices ni and ^2. The second DLM describes 
the transition from 712 to and the third DLM simulates the interface between media 
with refractive indices and n^. By inserting more DLMs before the rightmost one, it 
is possible to build any multilayer. For our simulation, we chose a periodic setup with 
alternating layers with reflective indices 712 and n^. 



i.e., if the stack ends on — )■ n^, and 



7^ 



N 



N-2\ 2 



1 



N-2 



if is even, 



(36) 



i.e., if the stack ends on n2 — )■ n^. 



5.2. Event-based simulation 

The method for simulating a plate consisting of two interfaces (section |4]) 
can be generalized to the case of a multilayer. In this case we concatenate N 
DLMs, each one simulating the behavior of a single interface. The schematic 
diagram of the DLM network is depicted in Fig. 15 A message, that is 



sent into the network can propagate back and forth through any DLM, until 
eventually it leaves the network through port of the first DLM or port 1 of 
the last DLM. This corresponds to multiple transmissions and reflections of 
a photon within the multilayer, until eventually it exits on the front or the 
back side. 
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Figure 16: Reflectivity 72. of a periodic multilayer as a function of the number of interfaces 
N. Closed circles: Event-based simulation data. Open circles: Wave theory (Eqs. (35) and 



(36l). The parameters were chosen to resemble an experiment with quarter wave films of 
zinc sulphide and cryolite at normal incidence [1]: ni = 1, n2 = 2.3, 713 = 1.35, 7x4 = 1.52, 
712/12 — "-3^3 — Ao/4, 9i — 0. The simulation was carried out with a ~ 0.998 and 1000000 
events per data point. The results of our event-based approach agree very well with the 
predictions of wave theory. 



5.3. Simulation results 

We study the case of a periodic multilayer. The sequence of refractive 
indices along the stack is ni, n2, ns, n2, na, . . ., n2, (^3), (figure [II)). 



Figure 16 shows the reflectivity (normalized intensity on output port 0) de- 
pending on the number of interfaces. Another analysis shows the reflectivity 



depending on the ratio ^2/7^3 (Fig. 17). 



5.4. Discussion 

We have shown that the approach introduced in section [2] can be gener- 
alized to multilayers by concatenating multiple DLMs. The basic building 
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Figure 17: Reflectivity TZ of a periodic multilayer depending on the ratio ri2/?^3 of the 
refractive indices for N = 5 (solid circles) and = 6 (open circles) interfaces. Data 
points are simulation results as obtained with a = 0.99 and 100000 events for each point. 
Solid lines are results derived from wave theory (Eqs. (351 and (36 1). Model parameters 
are: ni = 1, 712 = 2, = 1.52, 712/12 = n^h^ = Ao/4, 9i = 0, with varying. The 
event-by-event simulation gives the same results as the wave theoretical description (solid 
lines). 
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block is a DLM simulating a plane interface between two homogeneous dielec- 
tric media. A network of multiple DLMs can be used to form more complex 
structures. We studied the case of a periodic multilayer and found excellent 
agreement of the event-based simulation results with the wave theoretical 
description. 



6. Event-based simulation of a Mach-Zehnder interferometer 
6.1. Wave theory 

The Mach-Zehnder interferometer is a device that is sensitive to relative 



phase shifts [T]. The schematic setup is depicted in Fig. 18 At any time, 
there is at most a single photon traveling through the system. After passing 
the first beam splitter the photon gets an additional phase shift depending 
on the path that it follows. It passes the second beam splitter and is detected 
on one of the output ports. According to wave theory the amplitudes of the 
photons (&0)&i) in the output ports (A^2) and 1 (A^a) are given by 

6° ) = 2 ( 1 1 ) ( ( 1 ) ( «i ) ' 

with (ao,ai) being the amplitudes in the input ports and ipQ, ipi describing 
the additional phase rotations in the specific arm of the interferometer. 

For input on port only, i.e. (ao,ai) = (cos0o + i sin ^o, 0), we get the 
probability distribution 

|oo| = sm , \bi\ = cos . (38) 



6.2. Event-based simulation 

We simulate a Mach-Zehnder interferometer with our event-based ap- 
proach by using the same building blocks as shown in Fig. [18] we use two 
beam splitters, two phase shifters, one single-photon source (not shown) and 
two detectors (not shown). The beam splitters are built from more fun- 



damental single interfaces as described in section |4.4[ The communication 
between the components of the interferometer is mediated by single photons 
only. Each photon carries a message that is just its phase information. 

The simulation results for input on a single port only, i.e. (cos 0, sin 0) 



on port and no input on port 1, are shown in Fig. 19 The phase (j) is 



chosen randomly but fixed from an equal distribution in the range [0,360°]. 
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Figure 18: Schematic setup of a Mach-Zelmder interferometer, consisting of two beam 
splitters and two devices that perform a specific phase shift. A photon is sent on one 
of the input ports and leaves the beam splitter on one of its output ports. Depending 
on which path of the interferometer a photon takes, its phase is changed by either ipo or 
ipi. It then interacts with the second beam splitter and eventually, a photon is detected 
on one of the output ports of the interferometer. In the simulation we denote the photon 
numbers on each path with A^o, -^i; -^2, and N^. The counts Nq and Ni give the number of 
particles that travel from the first to the second beamsplitter via one of the two pathways 
(but never via the two pathways simultaneously) 
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Figure 19: Normalized intensities in Mach-Zehnder interferometer (see Fig. 18 ) depending 
on the phase shift tpo for various values of the phase shift ■01 ■ The data points are simulation 
results using 100000 events and a = 0.98. Sohd circles: No/{No + Ni) for V'l = 0; Open 
circles: N2/{N2 + N3) for tPi = 0; Solid squares: iV2/(A^2 + N3) for i/>i = 30°; Open 
squares: iV2/(^2 + ^^3) for tpi = 240°; Solid triangles: N3/{N2 + N3) for Vi = 0; Open 
triangles: Ny,/{N2 +-/V3) for V-'i = 300°; The solid lines are the predictions of wave theory. 



Each data point represents 100000 events (A^o + A^i = A^2 + A^s = 100000) 
and a is set to a = 0.98. The angle of rotation ipQ is varying in steps of 10° 
and the normahzed intensities for the various ports are determined from the 
event-based simulation. 

The agreement between the event-based simulation and the wave theo- 
retical description is very good if a is close to 1. This has been shown for 
event-based simulations with a DLM model that simulates a beam splitter 
as a whole [3],|ll|5], not as a collection of interfaces, as is done in the present 
work. Thus, we have shown that DLMs describing beam splitters can be 
built up from more fundamental building blocks, namely single interfaces, 
illustrating the modularity of our simulation method. 
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6.3. Discussion 

Our event-based simulation of a MZI shows that this optics experiment 
can be simulated by building a network of a basic building block, the DLM- 
based machine that simulates an interfaces between two homogeneous media. 
The simulation results for the basic building block and the more complicated 
networks such as a plane parallel plate, a beam splitter, a periodic stratified 
medium, and a MZI are all in excellent agreement with the corresponding 
results of Maxwell's theory, demonstrating that our event-based approach 
is modular. Crucial for the DLM-based simulation to yield the results of 
Maxwell's theory is that the parameter a which controls the dynamics of the 
DLM is close to one P H [5] . 

7. Conclusion 

We have shown that basic optical phenomena such as reflection, refraction 
and interference can be simulated by an event-based, particle-like approach. 
Our computational approach has the following features: 

• It yields the stationary solution of the Maxwell equations by simulating 
particle trajectories only. 

• Material objects are represented by DLM-based units placed on a bound- 
ary of these objects, which in practice involves some form of discretiza- 
tion. Apart from this discretization, all calculations are performed 
using Euclidean geometry. 

• Unlike wave equation solvers, it does not suffer from artifacts due to 
the unavoidable termination of the simulation volume P|: Particles 
that leave this volume can simply be removed from the simulation. 

• Unlike wave equation solvers which may consume substantial computa- 
tional resources (i.e. memory and CPU time) to simulate the propaga- 
tion of waves in free space, it calculates the motion of the corresponding 
particles in free space at almost no computational cost. 

• Modularity: Starting from the unit that simulates the behavior of a 
plane interface between two homogeneous media other optical compo- 
nents can be constructed by repeated use of the same unit. 
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We believe that the work presented in this paper may open a route to 
rigorously include the effects of interference in ray-tracing software. For 
this purpose, it is necessary to extend the DLM-based model for lossless 
dielectric materials to, say a Lorentz model for the response of material to 
the electromagnetic field |2]. This may be done by simple modifications of 
the DLM update rule, an extension that we leave for future research. 
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